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Abstract 

This article contains a review of an alternative theory of squeezing, 
based entirely on the wave function description of the squeezed states. 
Quantum field theoretic approach is used to describe the squeezing of 
the electromagnetic field in its most complete form that takes into ac- 
count temporal and spatial characteristics of the squeezed state. An 
analog of the Wigner function for the full electromagnetic field is in- 
troduced and expressed in terms of second order correlation functions. 
The field-theoretic approach enables one to study the propagation of 
the "squeezing wave" in space-time. A simple example of weak squeez- 
ing, that allows for all calculations to be done analytically, is discussed 
in detail. 

PACS: 42.50.Dv, 12.20-m, 03.65. Ca 



1 Introduction 

The standard description of squeezing of the electromagnetic field ^, ||, [|] 
employs the annihilation and creation operators associated with the mode 
decomposition of the field operators. Subsequently, only one or two modes 
are singled out for the analysis of the squeezing phenomena. This procedure 
gives a very effective tool to explain most squeezing experiments. However, it 
will not be adequate if one attempts to describe the propagation of squeezing 
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in space-time, when many modes are participating in the process and the 
state of the field exhibits strong correlations between different modes. In 
order to accommodate such situations, I shall introduce a description of 
squeezing, started in Refs. || ^, f?j and continued more recently ||, that 
employs formal methods of quantum field theory and enables one to produce 
compact and closed expressions for the space-time characteristics of the 
squeezing phenomena. Very convenient tools in this analysis will be the 
wave functional and the Wigner functional of the quantized electromagnetic 
field. These objects require a specification of the reference frame since they 
depend on the field 3- vectors B and D. At a first glance such a description 
may look noncovariant and the author of Ref. || has fallen into this trap. 
However, as I show in Sec. |8], despite the absence of tensorial indices, our 
formulation is fully relativistically covariant. 

Since the electromagnetic field may be viewed as one grand, multidimen- 
sional harmonic oscillator, many properties of squeezed states can be clearly 
seen in a simpler case of an ^-dimensional harmonic oscillator in quantum 
mechanics. Using this example, I introduce in Sec. | the necessary formal- 
ism and I describe in Sec. || special properties of the squeezed states related 
to uncertainty relations for positions and momenta. Next, I solve in Sec. || 
the time evolution problem for squeezed states and I introduce in Sec. | the 
description of squeezed states in terms of the Wigner function. 

In Sec. ^, the formalism developed for squeezed states of the quantum 
mechanical oscillator will be carried over to quantum electrodynamics. Even 
though the electromagnetic field corresponds to an infinite- dimensional os- 
cillator, there is a complete formal correspondence between these two cases. 
This correspondence enables one to write complete expressions for the rele- 
vant quantities in quantum electrodynamics in a closed form even though the 
rigorous mathematical foundation for this theory is often lacking. I address 
In Sec. [7| the problem of the propagation of squeezing of the electromag- 
netic field in space-time and I discuss in detail simple example of a small 
perturbation. Finally, I introduce in Sec. || the Wigner functional for the 
full quantized electromagnetic field — an analog of the quantum-mechanical 
Wigner function and I discuss its relativistic properties. 

There is no universally adopted terminology concerning the squeezed 
states and some authors extend the notion of squeezed states to almost any 
state that exhibits some kind of squeezing. Throughout this paper I shall 
use the term squeezed states in a narrow sense to denote the states whose 
wave function is a Gaussian. Such states, as discussed in Sec. ||, saturate 
various forms of the uncertainty relations for positions and momenta. They 
also can be obtained, as is discussed in the standard theory of squeezing, 
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from the vacuum state by the action of operators that are quadratic in the 
annihilation and creation operators in the exponent [||, ||, |j . 



2 Multidimensional harmonic oscillator 
in quantum mechanics 

Squeezed states are usually defined in terms of annihilation and creation 
operators, but the description in terms of the wave function in coordinate 
representation or in momentum representation, which I have adopted here, 
is very useful for a general analysis: 

• It offers a geometric picture of the squeezing phenomenon that is easy 
to grasp intuitively 

• It leads to compact formulas that exhibit the correlations between 
various oscillators 

• It enables one to give a simple general description of time evolution 

• It uses the smallest possible number of parameters to describe an ar- 
bitrary squeezed state 

Of course, the formalism based on annihilation and creation operators is 
often very convenient but for the uniformity of my exposition I shall restrict 
myself solely to the description in terms of wave functions. 

The Hamiltonian of an iV-dimensional harmonic oscillator will be taken 
in the form (summation convention for repeated indices will be used through- 
out) 



where both matrices g and u are symmetric and positive definite. One might 
also include the quadratic terms that mix x and p but I shall omit them since 
they are not needed in the analysis of the electromagnetic field. Changing 
the variables, one may transform g %3 to a multiple of the Kronecker delta. 



By an additional orthogonal transformation, one may bring the potential 
matrix u to a diagonal form. The eigenvalues itj of this matrix are multiples 





g 



m 
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of the squares of the characteristic frequencies fij of the oscillators, U{ = 

The ground state of the quantum-mechanical system with the Hamilto- 
nian (Q) can be obtained by simply plugging in a Gaussian wave function of 
the form 



V>o(r) = (DetO) 1 / 4 (m/7r^) Ar / 4 exp(-^x i fi ij x J ), 



(3) 



into the time-independent Schrodinger equation Hip = Eip and by solving 
the matrix equation 



m ftgft = u. 



(4) 



The solution for Q must be chosen as a symmetric and positive matrix. In 
our diagonal representation f2 has the form 



n 



/ fii 
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n n J 



(5) 



In addition to the ground state, there exist other Gaussian wave functions. 
The most general normalized Gaussian wave function of an iV-dimensional 
harmonic oscillator is characterized by two real vectors £ l and 7Tj and two 
real symmetric N x N matrices forming one complex matrix K^j = aij+ibij, 



ip(r) = (Deta) 1/4 (m/vr/i) JV / 4 exp(-^V i /2/i) 

TYl 

x exp[- — (x J - OKijix 3 - £ J ) + iinx % /h], 



(6) 



where the phase factor has been introduced to maintain a symmetry between 
the position @ and the momentum representation (p^). 

Gaussian wave functions are special: the Gaussian form of the initial 
wave function of a harmonic oscillator remains Gaussian during the whole 
time evolution. This property has been recognized right after the discovery 
of wave mechanics |jlC|]. In addition, as shown in Sec. ||, the Gaussian wave 
functions saturate the multidimensional generalization of the Schrodinger- 



Robertson uncertainty relations [11, |12|] and the entropic uncertainty rela- 



tions [13]. Gaussian distributions in probability theory are also exceptional: 
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all higher order correlation functions can be expressed in terms of the low- 
est correlation functions. This property is reflected in the properties of the 
Wigner function for the squeezed states, as shown in Sec. ||. 

The notion of a squeezed state always refers to a specified Hamiltonian. 
Given a Hamiltonian of the form (|l]) one can find the ground state but the 
inverse problem does not have a unique solution because for every symmetric 
and positive matrix there is a whole family of Hamiltonians of the form 
(|l|) such that the ground state is described by (Q). The only condition on 
the matrices g and u is that they satisfy the condition (Q). One can even 
find a Hamiltonian for which the ground state is described by a complex 
matrix K but if the imaginary part of K does not vanish the Hamiltonian 
must also contain the mixed term f i J x l pj. 

The choice of the ground state provides us with a standard against which 
one can measure the departures of the mean square deviations in positions 
and momenta from their normal (i.e. ground state) values. This enables 
one to define the notion of squeezing. In particular, it enables one to relate 
the scale in the momentum variables to that of the position variables. All 
Gaussian wave functions correspond to squeezed states, except the one de- 
scribing the ground state. When a is "larger" than f2, the state is squeezed 
in the position variables and when it is smaller, it is squeezed in the mo- 
mentum variables. Of course, since K is a matrix, the degree of squeezing 
may depend on the direction. 

The energy of the system separates into the center of mass part and the 
internal motion part 

E = E cm + E- mt , (7) 
E cm = ^^j + ^ee, (8) 

Emt = j Tr{mgba~ 1 b + mga + ua~ l /m}. (9) 

The energy of the internal motion tends to zero in the classical limit, when 
fi — > 0, because the Gaussian shrinks to zero and one is left with the mo- 
tion of a classical point-like particle. It is worthwhile to observe that the 
Gaussian solutions of the Schrodinger equation correspond to two indepen- 
dent autonomous mechanical systems with their Hamiltonians given by the 
expressions @ and @. The canonical variables for these two systems, as 
has been shown in Ref. [||, in addition to TVi), are the components of the 
matrices a~ 1 and b. 
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Prom the expression for Ei nt one may obtain the following formulas for 
the first and second variation of the internal energy around the ground state 

5E[ nt = — Tr{<5a — ua^daa' 1 /m}, (10) 

S 2 E int = ^TriSbn^Sb + Safl^Sa}, (11) 

where I have used g %3 = S v /m. The equilibrium condition derived from 
5E mt = (Sa is arbitrary) coincides with the equation (Q). The quadratic 
form defining the second variation is clearly positive confirming the fact that 
the energy of the ground state is the lowest one. 

The formula (0) can also be used to find the change in the energy of the 
ground state due to a change of the matrices g and u describing the system. 
Here one can see the advantage of keeping g and u arbitrary. 

% 

5E- mt = —Tr{mSg a + mg6a + Sua" /m — ua~ 5aa~ /m) 
h 

= — Tt{m5g a + 5u a _1 /m}. (12) 

Since for the ground state the variations 5a of the matrix a cancel out, this 
formula may be viewed as an analog of the Feynman-Hellmann theorem 

MM- 

Analogous results are obtained in momentum representation. The gen- 
eral normalized Gaussian wave function in momentum space is obtained by 
the Fourier transformation of the position space wave function (^) and it 
reads 

0(p) = (Deta)- 1/4 (vr/im)- iV/4 exp(^V i /2/i) 

x eM-^iPi-nKK^fipj-^-iCpi/h}. (13) 

Also the formula @ for the internal energy of the Gaussian state, despite 
its appearance, exhibits a symmetry under the interchange of positions and 
momenta. Denoting by c and d the real and the imaginary parts of the 
matrix K^ 1 , 

c = (a + ba-H)' 1 , d = -a" 1 b(a + ba- 1 b)- 1 , (14) 
one can rewrite ([?]) as 

Eiisk = —Tr{udc~ 1 d/m + uc/m + mgc^ 1 }. (15) 
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3 Saturation of uncertainty relations 
by squeezed states 



The special character of squeezed states comes forth clearly in the analysis of 
uncertainty relations between positions and momenta. In order to simplify 
the notation, I shall assume in this section that the expectation values (x 1 ) 
and {pi) are zero. In the simplest, one-dimensional case the Schrodinger- 
Robertson uncertainty relation |ll], |l2| has the form 

Ax 2 Ap 2 > h 2 /A + (xp + px) 2 /A. (16) 



This inequality is saturated by any Gaussian wave function. The equality 
expressing the saturation of uncertainty relations in one dimension can also 
be rewritten as a matrix equation 

QJQ = {h/4) 2 J, (17) 



where 



n-t \{xp + px) , 



and J is the simplectic metric 

The set of conditions that correspond to the saturation of the Schrodinger- 
Robertson inequality in the iV-dimensional case has the same general form 
( |l7| ) except that all matrix elements of the matrices Q and J become now 
N x N matrices 

Qn=( i r iff. iM m ^ P f ] ), (20) 

^ ^{piX 3 +X ] pi) {piPj) J 

and 

HA?)- 

In order to prove the formula (O) in the iV-dimensional case, one may use 
the following expressions for the position and momenta correlations func- 
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tions that can be obtained from the formula (g) by performing the corre- 
sponding Gaussian integrals 



(PiPj) 
(pix j ) = 



(x t x :l ) 
hm 



2m 



(a + ba b)ij, 



2 (5) + i(a~'b))), 



i(ba 



'3' 



(22) 
(23) 
(24) 
(25) 



These relations can easily be inverted to give K and its inverse in terms of 
correlations 



K ik = -i((xV)) 1 {x J p k )/m, 
(Kik)' 1 = im({p i p j ))- 1 (p j x k ). 



(26) 
(27) 

The validity of Eq. ( |l7j ) with Q and J given by the expressions ( 20|) and 
( ^l|) can be checked now by a direct calculation with the use of (^2 |25|). 

By a direct calculation one may also check that the entropic uncertainty 
relation fll3 | 

d N x\ 



fr)| 2 ln| 



A' 



P)| 2 ln| 



(p)| 2 >iY(lnvr + l) (28) 



is saturated by all squeezed states. Note that the saturation of the Schrodin- 
ger-Robertson and the entropic uncertainty relations is a characteristic fea- 
ture that singles out all squeezed states. The whole set of squeezed states 
can be characterized by these properties without any reference to a specific 
Hamiltonian or to annihilation and creation operators. 



4 Time evolution of squeezed states 

The squeezed states retain their form during time evolution. In order to 
satisfy the time-dependent Schrodinger equation, the vectors £, ir, and the 
matrix K must be taken as functions of time, i.e. 

i/)(r,t) = (Deta(t)) 1/4 (m/^) Ar / 4 exp(-if7r i /2n) (29) 

777 

x exp[-^(x»-e» (t)Kij(t))(^ -e(t)) + iwi(t)x*/h]. 
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Upon substituting this formula into the Schrodinger equation with the Hamil- 
tonian (||), one finds the classical equations of motion for the position £ and 
momentum ir of a harmonic oscillator 



£'(t) = s«>r,-(t), - r Mt) = -n ij e(t), (30) 



and the following nonlinear evolution equation for the matrix K 

^K(t) = -iK(t)mgK(t) + iu/m. (31) 

The vectors £ and it determine the center of mass position and the center of 
mass momentum of the Gaussian wave packet. The matrices a and b deter- 
mine the shape of the wave function and the distribution of the probability 
current of the internal "tumbling" motion, respectively. It follows from Eq. 
@ and from the relation dK~ x /dt = -K^dK/dtK' 1 that the inverse 
matrix K~ l satisfies the evolution equation in which the roles of g and u 
are interchanged, 

— K~ l (t) = -iK~ l (t) u/m K~ l (t) + img, (32) 
dt 

as one would expect from the formula (|l^) and the symmetry of the Hamil- 
tonian under the interchange of positions and momenta. 

The nonlinear equation ([H]) is a matrix Riccatti equation and it may 
be converted into a set of linear equations by the following substitution 
(g = l/m) 

K(t) = D~ 1 (t)N(t)n. (33) 



Eq. (31) is satisfied when two complex matrices N and D satisfy the set of 
equations 

^-D(t) = iN(t)n, ^N(t) = iD(t)Q. (34) 

til/ Ci(/ 

The solution of this linear set of equations 

N(t) = JV(0) cos fit + iD(0) sin nt, (35) 
D(t) = £>(0)cosfit + iiV(0)sinfit, (36) 

gives the following result for K(t) 

K(t) = (cosfit + ^(0)^^)" 1 (^(0)cosfit + ifisinfit). (37) 
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The order of matrices is here important since in general the initial value 
if (0) and fi do not commute. Actually, there are two equivalent forms of 
K(t) resulting from the identity 

(cos fit + iK(O)^^) -1 (if (0) cos fit + in sin fit) 

= (cos fit K(0) + ifisinfit)(cosfit + i^^^if(0))~ 1 . (38) 

Since both sides of this equality are related by the matrix transposition, the 
symmetry of the matrix K{t) is preserved during time evolution. For the 
ground state (and only for the ground state), K is constant in time, even 
though both matrices iV and D vary with time: N(t) = D(t) = e lf1< . From 



Eq. (37) one obtains the following formula for the time evolution of a 

a(t) = D-^ao-D 1 " -1 ^), (39) 

that shows explicitly that the positive definite character of a is preserved 
during the time evolution. 

The formula (37) describes completely the dynamics of squeezed states 
for an arbitrary harmonic oscillator. The time evolution of K may get quite 
complicated for a multidimensional harmonic oscillator when the matrices 
fi and if (0) do not commute so that they cannot be simultaneously diag- 
onalized. In the simple case, when the initial value of if represents just a 
small perturbation of the ground state, 

K(0) = Q + 5K(0), (40) 

the formula ( |37"| ) simplifies as follows 

K(t) = fi + e- int 5K(0)e~ int . (41) 

An alternative method of finding the time evolution will be given in the next 
section in terms of the Wigner function. 

5 Wigner function of squeezed states 

The Wigner function W sq (r, p) of a general squeezed state (||) can be easily 
determined from the defining formula [|l(| 

W(r,p) = 2~ N J d N riexp{iri-p/h)iP(r - rj/2)ip*(r + ri/2) (42) 
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by evaluating an iV-dimensional Gaussian integral. The prefactor in ( |42| ) 
has been chosen so that the Wigner function is a pure exponential, 

W sq (r,p) = exp[~(x - £)> + ba^b^x - 

x exp[-±(p ~ vUa-'fip ~ *)j -\{x- iy{ba~ l )i(p - vr),]. (43) 

This formula can be rewritten in a compact form in terms of the expectation 
values of the position and momentum operators 

^ sq (r,p) = exp[-^<(V " ~ (ft " <Pi»^) 2 >]- (44) 

Thus, the expectation values of the position and momenta and their quadratic 
forms determine the Wigner function. This was to be expected since ev- 
ery Gaussian distribution is determined by its second moments. For the 
squeezed states centered at the origin, this formula takes on the form 

2 

W sq (r,p) = exp[-~2 [x l (pipj)x J + pi{x l x 3 )pj - Pi(x l pj +p j x l )x 3 ]]. (45) 
n 

The time evolution of the Wigner function of a harmonic oscillator is deter- 
mined by the equation 

d t W(r, p, t) = -( Pi g ij V 3 - x L Uij d j )W{v, p, t), (46) 

where V and d denote the differentiation with respect to r and p, respec- 
tively. The solution of this equation can be written in the form 

W(r, p, t) = W(r(-t),p(-t),t = 0), (47) 

where (r(— t), p(— t)) is the time-reversed solution of the classical equations 
of motion. Owing to the invariance of the simplectic form x l pi — p{X l under 
all canonical transformations including the time translation, the formula 
(f45l) for the time-dependent Wigner function of a squeezed state can also be 
written in terms of the solutions of the Heisenberg equations of motion for 
the position and momentum operators 

VMr,p,i) = exp[--|<(y " <**(*)»&(*) " (Pi ~ <Pi(t)»^(t)) 2 >]- (48) 
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6 Quantized electromagnetic field 



Owing to the correspondence between a set of harmonic oscillators and 
the electromagnetic field, an analogous approach to the one presented here 
works, at least at the formal level, for the general squeezed states of the 
electromagnetic field. 

In order to use the analogy with quantum mechanics of harmonic os- 
cillators to its full extent, I shall employ the formulation of quantum elec- 
trodynamics in a representation which is the closest possible counterpart of 
the Schrodinger representation. Jackiw, who is also fond of this approach, 



called it in Ref. [17] the field theoretic Schrodinger representation. In this 
representation the state is described by a functional \& [A] of the vector 
potential. The vector potential operator A(r) and the magnetic induction 
operator B(r) acts on the state functional as a multiplication, 

A(r)# [A] = A(r)tf [A] , B(r)¥ [A] = (V x A(r))tf [A] , (49) 

while the electric displacement operator D(r) acts as the functional differ- 
entiation, 

D(r)*[A]=^^-^[A]. (50) 
The Hamiltonian in this representation takes on the form 

H = I fd 3 r [- f . - /. 2 . _ + -^(V x A(r)) 2 ]. (51) 



2 J L e(r)5A(r) 2 fi(r) 

This representation is the analog of the Schrodinger (position) representa- 
tion. There is also the analog of the momentum representation in which the 
electric displacement vector acts as a multiplication and the vector potential 
acts as the functional derivative. 

The most general Gaussian functional is determined by two vector func- 
tions A(r) and T>{r) and — the counterparts of £ and n — and a complex 
symmetric kernel /C(r, r') — the counterpart of the matrix K, 

*[A] = C&v[-^J<Prf<Pr'(A(r)-A(r)yK;{ryy(A(i f )-A(i'))] 
xexp[^- jd 3 r Z>(r)-A(r)], (52) 

In addition to being a function of two vector arguments, the kernel K, is also 
a 3 x 3 matrix since it acts on the vector indices of Ai. The vector function 
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A(r) and the kernel tC are just labels characterizing the state. To distinguish 
A from the argument A of the state functional, they are set in a different 
font. In order to secure gauge invariance, i.e., the invariance of ^ [A] under 
the gauge transformations A(r) — > A(r) + VA(r), the kernel X? J (r, r') must 
obey the conditions 

diK ij (r,r') = = K, ij (r,r')d'y (53) 

Thus, the kernel KL must be a double curl of some new kernel W, 

^i( r y) = e lkm d k W mn (r,r')d'^ ln . (54) 

This means that ^ [A] depends on A(r) only through B(r) = V x A(r). 

The total energy of the electromagnetic field, like in the case of the 
quantum-mechanical oscillator, is made of two parts: the classical energy of 
the coherent field and the quantum correction due to squeezing 

E ,u^m + m (55) 



2 J \ e(r) M (r) 
+ \ Jd'rjd'r' 6 ' 16 ^ r0 e imk d m a^(r, v')e 3nl d' n 

where a 13 and b lJ are the real and the imaginary part of K} 3 , 

V 3 (r, r') = a ij (r, r') + ib ij (r, r') (56) 

and the inverse of a 13 is to be taken in the subspace of transverse kernels 
satisfying the conditions (|53|), 

The analog of the equation (|j) that determines the ground state of the 
harmonic oscillator in quantum mechanics has now the form 

|d 3 rr fc (r',r,t)^/C^(r,r",t) = ^ J^'j^ d'^ 1 . (57) 

In the simplest case, when the medium is homogeneous, the kernel /C may 
depend only on the difference of coordinates and the solution of Eq. ( |57| ) 
can be easily found by the Fourier transformation. The resulting expression 
has the form 

V 3 (k) = fl e imk k m ^k n e 3nl . (58) 



13 



JC ij (r-r') 



1 1 _ i e 



I Jin 



(59) 



The coordinate representation of this expression is 

1 ft ' Akrnp, 

Therefore, the ground-state functional of the electromagnetic field in the 
homogeneous medium is 

1 IT /" o /".,,_, s 1 



[A] = C exp 



4vr 2 ?i 



- ld 3 rld 3 r'B(r) 



r r\2 



B(r') 



(60) 



This formula in the case of the vacuum state has been written down by 



Wheeler [18]. Here I would like to study its significance in the presence of a 
medium. The existence of the ground state in any static medium enables one 
to define the notion of squeezing with respect to this ground state. One can 
see from ( |60|) that the dielectrics decrease the fluctuations of the magnetic 
field, while the magnetics increase them. 

In the case of the electromagnetic field there also exists an analog of 
the quantum- mechanical momentum representation. In this representation 
the state of the electromagnetic field is given as a functional of D which is 
the canonically conjugate variable to the potential A. Owing to the known 
symmetry of Maxwell theory without sources, the formula for the ground 
state functional in the D representation can be obtained by the replacements 



D -! 

and it has the form 



B, B 



D, e 



e. 



C'exp 



(61) 
(62) 



where A is the vector potential for D, i.e. D = V x A. 

For an inhomogeneous medium, one may easily find first-order changes 
in K, by perturbation theory. The most interesting case is that of a change 
in the dielectric constant. The exact kernel for the ground state in a static 
medium satisfies the equation 



Urn / ,\ O ry ) 



(63) 



e(r) /i(r') 

from which the following equation is obtained by taking a variation due a 
small change in e(r) 

fd 3 r5V k (r',r)^-lC l i(r,r") + [d 3 r IC ik (r' ,r)^-5IC l ^r,r' f ) 
J e(r) J e(r) 

4z#e(r) 



d 3 rK. ik (r',r) 



e(rf 



-K, lj { 



r, r 



(64) 
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When the variation is taken around the vacuum state, the translational 
invariance of this state allows one to seek the change in fC in the form 

5)C ij (r',r") = Jd 3 rT ij (r' - r,r - r")8e(r). (65) 



Applying the Fourier transformation to Eq. fl64|), one finds the following 
Fourier transform of the kernel V 13 

r«(ki, k 2 ) = [n\ - - ks*i), (66) 

where c is the speed of light in the medium. In the coordinate representation, 
the expression for reads 

r«(r,.r 2 ) - 0»A - (67) 

This expression shows the long range effect in the squeezing produced by 
a dielectric: when one moves away from the dielectric, the function V falls 
only as 1/r 4 . 

With the use of Feynman-Hellmann theorem (^) one may directly arrive 
at the first order correction to the ground state energy Eq that is due to a 
small departure of the dielectric constant from its vacuum value eo, without 
ever calculating 5K. To this end, I shall use the following adaptation to the 
present case of the formula ( |T2] ) for the harmonic oscillator 



SE = ~ [d 3 r [cfr'^-SijK^ {r' -r,r-r') 
J e 



4 



1 f 3 5e(r) n f d 3 k 



= -iF r ^ L2 Jwr nc ^ (68) 

The integral over k represents the energy of all photons that are affected 
by the change in the dielectric constant and the factor of 2 comes from 
two polarization states of the photon. Obviously, the integration does not 
extend to infinity because all dielectrics become transparent (behave like the 
vacuum) for sufficiently high photon frequency and that fact will introduce 
a cutoff. Therefore, a more appropriate way to write Eq. ( |68|) is 

5Eo = -\Jd*r2j£L^HcUr)\K (69) 

where e and c are taken as slowly varying functions of the wave vector k. 
Since 5c = 5(l/y/ejT) = — (Se/2e)c, the expression (|69| ) can be recognized as 
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the variation of the zero point energy evaluated in the local approximation 
when the variation of e in space is very slow over the distance of a wave 
length 1/k, 

5E = 5(- jd 3 r2j-^ ftc k (r)|k|). (70) 

First order perturbation theory, described here as a simple illustration 
of the global approach to squeezing, is clearly not sufficient to give the 
Casimir force. This force can be understood as coming from the change in 
the energy of the ground state but it requires the calculation of the energy 
at least up to second order. In the first order, one obtains only the volume 
effect. In particular, the change of the relative positions of several pieces of 
the dielectric would not lead to a change in the energy and therefore would 
not give rise to an interaction force. 



7 Propagation of squeezing 

In order to satisfy the time-dependent Schrodinger equation 

ihd t V [A\t] = [A\t] , (71) 

the functions A and T> and the kernel KL must depend on time, except for 
the ground state. Their time evolution can be derived by substituting the 
wave functional 

*[A|i] 

= C exph^r Jd 3 r Jd\' (A(r) - A(r, t)) -K(r, r', i) • (A(r') - A(r', t))] 
+ exp[^ Jd 3 rV(r,t)-A(r)] (72) 



into (|7l| ) and then by comparing the terms of the same form on both sides of 
this equation. In this way one obtains the Maxwell equations for B = V x A 
and T> 

^(r,() = -Vx^, (73) 
P(r,t) = Vx^M, (74) 
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and the following nonlinear integro-differential equation for K, 

d t Vi(r',r",t) (75) 
= -i jd\ K ik {v', v, t) ^JC'J (r, r", t) + ie imh & m " ^ d'^ nl . 

Propagation of squeezing, in principle, is fully described by the time 
evolution equation (|37| ) of the kernel JC but it is rather difficult to draw 
physical conclusions from this equation in the general case. However, for 
small perturbations, when the departure from the vacuum situation is given 
by (|65|) and ( p7| ) , one can find an explicit expression for the time dependence 
of /C. To this end, I shall use the first order formula (fill) obtained for 
the harmonic oscillator. In the present case, this formula takes on a very 
simple form when K, is taken as a function of wave vectors because in this 
representation the frequency O is diagonal. The application of the operations 
exp(— iClt) in the formula (41) amounts simply to the multiplications by 
exp(— ic|k|t). This leads to the following formula for SK(t) 

^'(ri,r 2j t) 

= J d \ y'^^e ikl < ri - r )e ik ^ r - r2 )5e(r) e --l kl l i e- ic l k2 l t f^(k 1 ,k 2 ) 
ftfVr^'(n -r,r-r 2 ,i)<5e(r), (76) 



The function T(ri — r, r — r 2 , t) can be easily evaluated explicitly since after 
the substitution of the expression ( |66| ) all integrations over the wave vectors 
are elementary. The resulting formula reads 

r^(n-r,r-r 2 ,t) 

= c(8 ik A - &&) x 8u,G{r x - r, r - r 2 , t)(^ A - d 1 8 j h, (77) 

where 

G(nr.t) - [ Al# \ta e ^,J^!l (78 ) 

G(ri ' r2 ' t} - ]j^fWf e e p#|nTP^) (78) 

1 f°° „ f°° „ e- icklt e~ ick2t snxVkxn) sin(fc 2 r 2 ) 



(2ir 2 ) 2 Jo Jo h + k 2 n r 2 

The integration over k x and k 2 can be easily performed after the change of 
variables: k = k\ + k 2 , q = k\ — k 2 

G(ri,r 2 ,t) = -— = / — e lKa \ dqsm——r 1 sm——r 2 



(27r) 3 7rrir 2 J Q k J-k 2 2 
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dk 



-ikct 



(2vr) 3 7r(r2 - rj) 
1 

(27r)3(r?-rg) 
Kr2~\ct\) 9( n 



sin kr2 sin fori 



''2 



'1 



(79) 



ct 



n 



H In 

7rr 2 



r 2 - ct 



r 2 + ct 



i 

7rri 



In 



ri — ct 



ri + ct 



where 0(r — ct) is the step function. The function (fn|) describes fully the 
propagation of squeezing in the vacuum when the departure from the vac- 
uum state is small (cf. Eq. (|4l|)). Unfortunately, this approximation is not 
valid near the light cone in the variables r and t because of the singularities 
in Eq. (|79|). A better, nonperturbative approach is needed to determine 
the behavior of squeezing on the light cone. What does follow, however, 
from these considerations is that squeezing has very interesting propagation 
properties. In my simple example, the initial value of the squeezing kernel 
K, has been influenced by a change in the value of e that took place in some 
region of space. This influence extends throughout space with the power- 
law falloff when moving away from the region where e has been changed, as 
determined by the formula (|67|). Then, at t = 0, the vacuum value of e is 
restored and the electromagnetic field begins its relaxation to the true vac- 
uum state. The real part of K, relaxes to zero with the speed of light as seen 
from the presence of the step functions, while the imaginary part relaxes to 
zero much more slowly as determined by the logarithmic functions. 



8 Wigner functional for the electromagnetic 
field and relativist ic invar iance 

There is a tendency among some physicists to limit the notion of relativistic 
invariance to theories that involve only tensorial (or spinorial) equations. 
The case of the Wigner functional and the number of photons, that is ex- 
plored in this section, shows clearly that this limitation is not justified and 
that it restricts unnecessarily the freedom of choosing the right formula- 
tion. All that is really required of a relativistically invariant theory is that 
the physical laws are the same in all inertial reference frames. This implies 
that the equations used by all observers related by a Poincare transforma- 
tion have the same form. The exact nature of the transformation laws for 
physically relevant quantities is of no importance if they only guarantee the 
correct final result. There is no simple way to write the expression for the 
Wigner functional for the vacuum state in terms of tensor quantities with 
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properly summed up indices even though this functional is an absolute in- 
variant. It has the same value and the same form for every observer and it 
is even invariant |23[| under conformal transformations. 

The description of squeezed states of the electromagnetic field in terms of 
the Wigner function can be obtained by analogy with quantum mechanics. 
As a starting point of this generalization I shall take the formula (fij|). The 
generalization consists of replacing the position r by the vector potential 
and the momentum p by (minus) the electric displacement vector. In this 
way, one obtains the following functional of A and D describing the most 
general squeezed state of the electromagnetic field 



WWA,D] 

= eM ~k 


jd 3 r 


fd 3 r' (A(r) 


-(A(r)))-(D(r)D(r'))-(A(r') 


(80) 
- <A(r')»] 


xexp[-A 


fd 3 rj 


fd 3 r' (D(r) 


-(D(r)))-(A(r)A(r'))-(D(r') 


-<D(r')»] 


xexp[-^ 


fd 3 rj 


fd 3 r' (A(r) 


-(A(r)))-(D(r)A(r'))-(D(r') 


- (D(r')))] 


xex P [-A 


jd 3 rj 


fd 3 r' (D(r) 


-(D(r))).(A(r)D(r'))-(A(r') 


-<A(r')»]. 



Owing to gauge invariance, the Wigner functional depends on the vector 
potential A only through the magnetic induction vector B and I shall take 
this explicitly into account by treating W as a functional of B and D. It 
should be stressed that the Wigner functional for the electromagnetic field 
is quite different from the Wigner function for the photon wave function 



discussed in Refs. pO| , 21]. The Wigner functional describes the quantum 
statistical properties of the full electromagnetic field and not just of a single 
photon. 

In the simplest case of the vacuum state of the electromagnetic field the 
expectation values of field operators can easily be evaluated and the resulting 
formula for the Wigner functional can be written in a form resembling the 
functionals representing the vacuum state vectors ( |60|) and (|62| ) 



exp(-2JV[B,D]), (81) 
where 

iV[B,D] (82) 
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The physical significance of the functional N [B, D] has been discovered by 
Zeldovich long time ago. It is the average number of photons present 
in the electromagnetic field described by the vectors B and D. This inter- 
pretation can be justified by noting that the number of photons is related 
to the energy of the field by the factor l/?ic|k| and the division by |k| is 
represented in the position space by a nonlocal kernel 

d 3 ke iMr~r')_L = 1 (83) 

|k| 47r 2 |r-r'| 2 

The functional (]82] ) that made its appearance in the exponent of the Wigner 
functional for the vacuum state has also been used as a natural norm for 
the electromagnetic field and it is intimately connected with the concept of 



the photon wave function [21, 22|. Despite its "nonrelativistic" appearance 



iV[B,D] is invariant not only under all Poincare transformations but also, 



as shown by Gross [ 23 1 , under the full conformal group. The formula (|82| ) 
teaches us an important lesson on the subject of relativistic invariance that 
appearances are sometimes misleading. 

A general Wigner functional (|8TJ) is, of course, not invariant under the 
Poincare transformations, but it does have well defined transformation prop- 
erties that follow from its physical meaning. As a (pseudo) probability dis- 
tribution in the phase space, the Wigner function and the Wigner functional 
must transform as scalars, because the probability is a scalar and the phase- 
space volume is an invariant. This requirement is even more obvious for 
the squeezed states, when the Wigner functional is a true probability dis- 
tribution. Thus, under a Poincare transformation, the Wigner functional 
transforms as follows 



W' [B', D'] = W [B, D] . (84) 

The transformed functional evaluated for the transformed values of its ar- 
guments is equal to the original functional for the original values of the 
arguments. The relativistic invariance of the theory means that there exists 
a set of generators of the Poincare transformations that obey the commu- 
tation relations characteristic of the Poincare group. These generators are 
identified with the corresponding physical quantities: energy, momentum, 
angular momentum, and the center of mass. In our case these generators 
are given by the standard expressions for the electromagnetic field and they 
act on the Wigner functional through the Poisson bracket relations. For ex- 
ample, the energy of the electromagnetic field generates the time translation 
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and that leads to the time evolution equation for the Wigner functional 



d t W [B, D|t] = {\ jd\ + [B, D|t]} (85) 



that is the counterpart of (t46|). The complete analysis of relativistic in- 
variance along these lines has been carried out in detail in Refs. [24, plf| . 



The second reference is particularly relevant since it deals with the classical 
statistical theory of the electromagnetic field and all the formal properties 
of the Wigner functional are the same as those of the classical distribution 
function. 
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